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High contrast matter-wave interferometry is essential in various fundamental quantum mechani¬ 
cal experiments as well as for technical applications. Thereby, contrast and sensitivity are typically 
reduced by decoherence and dephasing effects. While decoherence accounts for a general loss of 
quantum information in a system due to entanglement with the environment, dephasing is due to 
collective time-dependent external phase shifts, which can be related to temperature drifts, mechan¬ 
ical vibrations or electromagnetic oscillations. In contrast to decoherence, dephasing can in principle 
be reversed. Here, we demonstrate in experiment and theory a method for the analysis and reduction 
of the influence of dephasing noise and perturbations consisting of several external frequencies in 
an electron interferometer. This technique uses the high spatial and temporal resolution of a delay 
line detector to reveal and remove dephasing perturbations by second order correlation analysis. It 
allows matter-wave experiments under perturbingly lab conditions and can be applied in principle 
to electron, atom, ion, neutron and molecule interferometers. 


I. INTRODUCTION 

Matter-wave interferometers accomplish an outstand¬ 
ing sensitivity for high precision measurements. This has 
been shown in various fields, e.g. by the measurement 
of inertial forces mig, the gravitational acceleration [3] 
and atomic or molecular polarizabilities HE]. As a con¬ 
sequence of the high sensitivity, interferometers are also 
susceptible to external perturbations and noise that de¬ 
phase the matter-wave. Typical dephasing sources are 
temperature drifts, mechanical vibrations or electromag¬ 
netic oscillations. These influences result in a ’’washing- 
out” of the interference fringes. It generally reduces 
the interference contrast [6] and is especially an obstacle 
in precision phase measurements such as in Aharonov- 
Bohm physics m m. Dephasing also needs to be pre¬ 
vented in decoherence studies where the transition from 
quantum to classical behaviour is analyzed through the 
contrast loss in matter-wave interferometers [Ml]. De¬ 
coherence destroys the phase relation between two beam 
paths by interaction with the environment, whereas de¬ 
phasing is due to a collective time-dependent phase shift 
of the wave functions. Both reduce the contrast and need 
to be distinguished to study experimentally the theory of 
decoherence [12]. The knowledge of decoherence mech¬ 
anisms is of importance to achieve long coherence times 
in a variety of hybrid quantum devices in the community 
to date dMH]. 

Gertainly, isolating, damping or shielding of interfer¬ 
ometers reduce the dephasing influences. However, in a 
perturbingly environment this often cannot be achieved 
efficiently in a broad frequency range. In addition, for 
some applications such as mobile interferometric devices 
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m, the perturbation sources and therefore the shielding 
requirements are changing between different experimen¬ 
tal runs. 

In this article we present a general method to reveal 
dephasing and reduce the influence of external multi¬ 
frequency perturbations and noise by second order corre¬ 
lation analysis of an interference pattern. The technique 
can in principle be applied in all matter-wave interferom¬ 
eters, where the detector has a high spatial and temporal 
resolution. Such detectors are available for electrons m, 
ions [20] , neutrons [21] and neutral atoms [22] . Perturba¬ 
tion frequencies below the mean particle count rate can 
then be revealed after data acquisition is completed. 

Recently, we demonstrated [23] that the influence of a 
single perturbation frequency can be isolated, removed 
and the complete spatial interference pattern be recov¬ 
ered. All parameters such as the contrast, pattern pe¬ 
riodicity, perturbation frequency and amplitude can be 
extracted from the perturbed interference data. This is 
useful to eliminate a distinct electromagnetic oscillation 
such as 50 Hz noise from the electric network or a vibra¬ 
tion from a rotating vacuum pump. 

In this article we enhance this method towards several 
frequencies with different amplitudes, that dephase the 
interference pattern at once. We provide a detailed theo¬ 
retical description of the multi-frequency case, for which 
we calculate an exact solution for the second order corre¬ 
lation function. While the theory holds for an arbitrary 
number of discrete perturbation frequencies, it requires 
increasing computational power. However, an approxi¬ 
mate solution can be found, which is used to extend the 
theory to arbitrary noise spectra with a high-frequency 
cutoff. Such noise is often found in typical lab situations, 
where several low frequency noise sources disturb the in¬ 
terference pattern at the same time. Here, the relative 
phases of the frequency components do not influence the 
correlation function and the approximate solution is ap¬ 
plicable. 
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We demonstrate our method by the application of an 
electron matter-wave interferometer in combination with 
a delay-line detector. Using an external magnetic field, 
we introduce artificial multi-frequency noise, which we 
analyze on the basis of the ” washed-out” interference pat¬ 
tern. In detail, we demonstrate two- and three-frequency 
perturbations, including the special case, where the fre¬ 
quencies are multiples of each other. The influence of 
broad-band noise is investigated by using a perturbation 
spectrum with a cutoff frequency of 150 Hz. For all sce¬ 
narios, the measured particle correlations are compared 
to our theory. 


II. THEORY 

Recently, we demonstrated a full correlation analysis 
for a single frequency perturbation of a matter-wave in¬ 
terference pattern in an interferometer [23]. In the fol¬ 
lowing we discuss the general case and extend the corre¬ 
lation analysis to multiple perturbation frequencies and 
broad-band noise. 

Similar as before, we start with the probability distri¬ 
bution f{y,t) of particle impacts at the detector 

f{y,t) = foil + K cos {ky + (p(t))) . (1) 


yielding 

5^^) (w, t) = 1 + ^K^ + e*''"A_] , (5) 

with 

N 

j=l nj,mj 

_ 0*[%’f)] f'j\ 

Th j j TTl/ j ^ \ / 

where Jn denote the Bessel functions of first kind and the 
sum in Eq. © has to be taken over all integer-duplets 
{nj^rrij} G Z, jf = 1... N, for which 

+mj)ujj =0 (8) 
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is fulfilled. Although this makes in principle for an in¬ 
finite number of addends, the contribution of integer- 
duplets {nj^rrij} > (j)j to A± is suppressed due to the 
strong decay of the Bessel functions. 

For multiple perturbation frequencies ujj and broad¬ 
band noise spectra, the constraint from Eq. © requires 
usually rrij = —rij for all j = 1.. .N. With Xnj,-nj — 
(—, this results in a correlation function 


Here K denotes the contrast, k = 27 t/X the wave num¬ 
ber of the unperturbed interference pattern with spatial 
periodicity A, and /o assures normalization. The pertur¬ 
bation is mediated via the time dependent phase shift 
which we assume as superposition of N discrete 
harmonic frequencies 


N 

(fit) = ^ ipj COS {ujjt + . (2) 

j=i 


Thereby ipj and (j)j denote the peak phase deviation and 
the perturbation phase at the frequency component ujj. 
These perturbations will typically lead to a ” wash-out” 
of the time-averaged interference pattern at the detector. 

To regain information about the unperturbed interfer¬ 
ence pattern, the spatial and temporal information of the 
particles arriving at the detector is recorded and evalu¬ 
ated by correlation analysis. Using Eq. Q, the second- 
order correlation function reads 


g^‘^\u,T) 


<^fiy + u,t + T)fiy, t) 

^ fiy + u,t + t) fiy, t) '>y,t 


, ( 3 ) 


with ^ denoting the average over position and 

time 


I pT rY /2 

« fiy, t) »,,= — fiy, t) dydt . 

( 4 ) 

In the limit of large acquisition times T ^ 2^juj and 
lengths Y ^ l//c, Eq. © can be solved analytically 


= 1A (r) cos (ku) , (9) 


with 


N oo 


Air) = E . (10) 


jf=1 nj = — oo 


JnA^jf cos(njujjr) 


The correlation function for multiple dephasing oscil¬ 
lations will thus show a periodic modulation in the spa¬ 
tial distance u between two detection events. It has the 
same periodicity as the undisturbed interference pattern. 
Therefore, it can proof matter-wave interference in a situ¬ 
ation where the periodic pattern would be ” washed-out” 
after signal integration. The amplitude A{r) of this mod¬ 
ulation, however, depends on the correlation time r, and 
is given by the specific perturbation spectrum. The func¬ 
tion A{r) includes a superposition of sidebands at dis¬ 
crete frequencies rijUJj {rij G Z). Their strengths are de¬ 
fined by the peak phase deviation ipj via the Bessel func¬ 
tions Jn^{(pj). Maximal modulation amplitude of O.bK^ 
is achieved at correlation time r = 0. 

While for N = 1, the results from Eq. ^ and (10), 
coincide with our previous result found for single fre¬ 
quency perturbations [23|, special care has to be taken 
for small numbers of perturbation frequencies, which are 
multiples of each others. In this case the constraint from 
Eq. (Isl) might be fulfilled for integers mj ^ lead¬ 


ing to additional terms in Eq. (10). This can be easily 


seen in the two-frequency case, with uj 2 = 2co’i, where 
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Eq. (|8| is not only solved for mi /2 = —nij 2 ^ but for 
all muTtiplets with ni + mi = —2 (n 2 +m 2 ), such as 
{ni,mi,n 2 ,m 2 } = {2,2,-1,-1}. Following Eq. @-0, 
additional terms appear in the correlation function, with 
introducing a phase shift to the cor¬ 
relation pattern, which depends on the correlation time 
and the perturbation phases (pj. The correlation function 
can then not be described by the simple form of Eq. 
but has to be calculated via Eq. (§-0. 

Similar effects may arise for larger numbers of pertur¬ 
bation frequencies, but the contribution of these parts 
become more and more negligible with respect to the 
main correlation component described by Eq. 

Examples for all of the above scenarios will be given in 
section lYl 


III. EXPERIMENT 


To demonstrate our method we use an electron biprism 
interferometer and analyze interference patterns that 
were subject to multi-frequency perturbations. The ex¬ 
perimental setup is equal to the one used for the single 
frequency correlation analysis and is described elsewhere 
[231427] . A coherent beam of electrons is field emitted 
from a single atom tip cathode [28l [29] and guided to¬ 
wards an electrostatic biprism wire. Being on a positive 
potential it separates and recombines the matter-wave 
leading to an interference pattern [30|. After magnifica¬ 
tion the interferogram is amplified by two multichannel 
plates and detected by a delay line detector with high 
spatial and temporal resolution [20]. The whole setup is 
in an UHV-chamber at a pressure of 1.4 x 10“^^mbar. 
The beam path inside the chamber is shielded by a mu¬ 
metal tube. 

For the course of this paper, we artificially disturb the 
pattern by several superposed oscillating magnetic fields. 
They are generated by two external magnetic coils in 
Helmholtz configuration positioned around the vacuum 
chamber. The frequencies are introduced by a frequency 
generator. The oscillating magnetic field shifts the in¬ 
terference pattern perpendicular to the fringes and the 
integrated pattern gets ”washed-out”. 

The extraction of the second order correlation function 
from the detector signal data was performed according 
to the same procedure as described in [23]. From the 
detector we get the position yi and the time ti of all 
particle events i = 1...M. We then determine the number 
Mu,t of particle pairs (i^j) with (yi—yj) G [u^ u-{-Au] and 
{ti — tj) G [r, r + Ar], yielding the correlation function 


g<-^\u,T) = 


TY 




M‘^ArAu 




( 11 ) 


Here, Ar and Au denote the discretization step size in 
correlation time and position, respectively. TYjM^ is a 
normalization factor and [{1—t/T){1 — \u\/Y)]~^ corrects 
Mu^t for the finite acquisition time T and length Y. 



FIG. 1: (Color online) Frequency dependent attenuation of 
the mu-metal shield around the beam path. The data was 
obtained by applying single frequency perturbations between 
0 and 500 Hz to the interferogram. Each point corresponds 
to a correlation analysis of a disturbed interference pattern 
where the extracted perturbation amplitude has been set in 
relation to the phase shift at constant magnetic field (0 Hz). 
The attenuation curve has been fitted by an appropriate func¬ 
tion. 


Before applying multi-frequency noise, the mu-metal 
shielding has to be characterized, as it causes a frequency 
dependent damping of the externally applied fields. This 
is done by applying a single frequency perturbation via 
the external coils with the perturbation frequency step¬ 
wise increased from 0 to 500 Hz. For each measurement 
the same oscillation amplitude was applied to the coils 
and an electron interferogram was recorded. The peak 
phase deviation of the perturbation in the individual in¬ 
terference pattern was extracted by fitting the correla¬ 
tion theory from Eq. (HI)-®, for the single frequency 
case {N = 1), to the measured particle correlations. The 
resulting phase deviations (^i{uj) were compared to the 
phase shift for constant magnetic field (/:?i(0), which can 
be extracted from the 0 Hz measurement and is uninflu¬ 
enced from the mu-metal shield. The frequency depen¬ 
dent attenuation of the mu-metal shield is then given by 
1 — Lpi{uj)/ipi{()) and shown in Fig. [^together with an 
appropriate fit function. 


IV. RESULTS 

For the experimental verification of our multi¬ 
frequency correlation theory we investigate different per¬ 
turbation scenarios. Starting with two and three pertur¬ 
bation frequencies, we also analyze a full perturbation 
noise spectrum. For all measurements an interference 
pattern with 7 x 10^ counts was recorded. 

In the first measurement the interference pattern was 
perturbed by two superposed magnetic oscillations at 
uji = 27r X GHz and ci ;2 = 27r x 40Hz. Using eq. 11, 
we then extracted the second order correlation function 
^expt('^ 5 ^) these data, which is shown in the up¬ 

per panel of Fig. [^ The theoretical correlation function 
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FIG. 2: (Color online) top: Experimental correlation function 
5'expt with two perturbation frequencies of = 27r x 6 Hz 
and 012 = 27r X 40 Hz. bottom: Theoretical fit according to 
Eq. 1^-0 with N — 2. 


5'thLr('^’ r), in the lower panel of Fig. 2, was fitted to the 
experimental data according to Eq. (5)- 0 with N = 2 
and shows good agreement with the experiment. From 
the theoretical fit we extracted the contrast K = 62.9 %, 
the pattern periodicity A = 2.08 mm and the peak phase 
deviations cpi = 1.34 tt and (p 2 = 0.93 tt. The phases 
01 and 02 have no influence at the resulting theoretical 
function in Fig. because the constraint from Eq. 0 
is only fulfilled for rrij = —rij and therefore Eq. Q be¬ 
comes Xnj.-rij ~ (—l)’^'^ - In this case we also could have 
used Eq. 0 -(prl to fit the experimental data, yielding 
the same result. 

With the obtained parameters from the theoretical fit 
it is possible to recalculate the interference pattern from 
the perturbed one in Fig. |^a). To get the unperturbed 
interference pattern we have to calculate the new posi¬ 
tions Unew according to the extracted perturbation fre¬ 
quencies ujj and peak phase deviations cpj 


X ^ 

Vnew = y - 

i=i 

with the particle positions of the perturbed pattern 
pattern periodicity A and the phases of the perturbations 
0j. As the phases could not be obtained from the corre¬ 
lation analysis, they have to be varied to find the maxi¬ 
mum contrast of the interference pattern. For our mea¬ 
surement we found 0i = — 0.337r and 02 = — 0.637r. The 
resulting back-calculated interference pattern is shown in 
Fig. ^b). In the dephased pattern only a small contrast 
of 6.4% can be adumbrated. After recalculating the new 
particle positions the interference pattern is clearly visi¬ 
ble and the contrast is significantly higher and amounts 
62.4%. 


(a) (b) 
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FIG. 3: (Golor online) (a) Electron interference pattern per¬ 
turbed by two frequencies cji = 27r x 6 Hz and UJ 2 = 27vx 40 Hz. 
The remaining contrast is A = 6.4%. (b) Recalculated in¬ 
terference pattern where the dephasing could be corrected 
by the presented correlation analysis together with Eq. (12) 
for N — 2. The contrast of the reconstructed pattern is 
A = 62.4%. 



EIG. 4: (Golor online) top: Electron interference measure¬ 
ment with two perturbation frequencies, that are multiples of 
each other, cji = 27r x 20 Hz and UJ 2 — 271 x 40 Hz. bottom: 
The exact solution from Eq. 0-0 with N — 2 was fitted to 
the experimental correlation data, showing good agreement. 


In a second experiment the correlation function was 
studied for the case of two perturbation frequencies that 
are multiples of each other. With uji = 27r x 20 Hz 
and co ’2 = 27r X 40 Hz the resulting correlation pattern 
is shown in Fig. As discussed in the theory section, 
the constraint from Eq. 0 is now not only fulfilled for 
mi /2 = —^ 1 / 2 - Thus additional terms appear in the cor¬ 
relation function, that result in a phase shift of the cor¬ 
relation pattern depending on the correlation time r and 
the perturbation phases 0j. This phase shift is clearly 
seen in the correlation pattern of Fig. In this case 
the full theory from Eq. 0-0 is required to fit the ex¬ 
perimental data, as the correlation time dependent phase 
shifts are not included in the approximate solution from 
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FIG. 5: (Color online) top: Measurement with three pertur¬ 
bation frequencies of oji = 27r x GHz, 002 = x 23Hz and 
6^3 = 27 r X 40 Hz. bottom: Fit from Eq. ([^-Q with N = 3 
to the experimental correlation data. 


eq. J^-(IO). The fit is shown in the lower panel of 
Fig. I^and shows good agreement with the experiment. 
From the fit we extract the contrast K = 62.9 %, the pat¬ 
tern periodicity A = 2.08 mm, the peak phase deviations 
(fi = 0.50 TT and (f 2 = 0-52 tt and the phases (pi = —0.21 tt 
and <p2 = 0.24 TT. 

Also the case of three superposed perturbation fre¬ 
quencies was studied and compared to our theoretical 
approach. The outcome of a corresponding electron in¬ 
terference experiment is shown in Fig.j^at the top. Here, 
the frequencies uJi = 27r x GHz, co ’2 = 27r x 23Hz and 
Co’s = 27r X 40 Hz have been applied. For fitting these 
data we use the exact solution from Eq. <§- 0 , for three 
frequencies, because the constraint from Eq. 0 is ful¬ 
filled for many multiplets with rrij ^ As 

for the case of two frequencies, that are multiples of each 
other (Fig. 1^, this leads to additional terms in the cor¬ 
relation function. The result of the fit is shown in the 
lower panel of Fig. yielding the contrast K = 61.3%, 
the pattern periodicity A = 2.08 mm, the peak phase de¬ 
viations ifi = 0.76 TT, ip 2 = 1.01 TT and = 0.52 tt and 
the phases pi = 0.01 tt, p 2 = 1.01 tt and = 0.01 tt. 
With the exact solution a good agreement between the¬ 
ory and experiment can be achieved. 

In the last experiment we investigate a perturbation 
consisting of a broad frequency spectrum. This should 
simulate typical noise from a lab environment originat¬ 
ing from electric or magnetic devices or mechanical vibra¬ 
tions. We introduced a perturbation noise spectrum from 
OHz up to a cutoff frequency at 150 Hz to the magnetic 
coils. A part of the time signal of the applied perturba¬ 
tion spectrum was measured while integrating the signal 
at the detector. The Fourier spectrum of this time signal 
(see inset of Fig. was attenuated according to the fit 
function shown in Fig. [^and the amplitudes of the fre- 


0.8 



UJjl2TT [Hz] 

FIG. 6: (Golor online) Perturbation spectrum at the position 
of the electrons which was used to calculate the theoretical 
correlation function in Fig. according to Eq. 

The spectrum originates from the externally applied noise 
spectrum (see inset) and the frequency dependent damping 
of the mu-metal shield (see Eig. [^. 



EIG. 7: (Golor online) Gorrelation analysis of an interfer¬ 
ence pattern that was perturbed by a broad-band noise spec¬ 
trum with a cutoff frequency at 150 Hz. This simulates 
a perturbingly lab environment, top: correlation function 
^expt('^ 5 '^) obtained from the measurement, bottom: The¬ 
oretical fit according to Eq. 


quency components at the positions of the particles were 
calculated. The result is shown in Fig. These values 
were used to calculate the theoretical correlation func¬ 
tion according to Eq. In this case it is possible 

to use the approximate solution for the correlation func¬ 
tion, because the perturbation consists of a large number 
of frequencies and uncorrelated phases. Fig. shows 
the measured data together with the theoretical fit, from 
which we deduce the contrast of the interference pattern 
K = 63.4% and the pattern periodicity A = 2.16 mm. 
Again the experimental and theoretical correlation func¬ 
tions agree well. 
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V. CONCLUSION 

Matter-wave interference experiments appear in vari¬ 
ous fields of experimental quantum physics for atoms m, 
molecules [32] , neutrons [33| , electrons [25] , and ions [26] . 
Most of these approaches have a spatial resolution for the 
interference pattern, but not a temporal one. Therefore, 
these kind of setups are susceptible to dephasing effects. 
Oscillations from electromagnetic sources, temperature 
drifts or mechanical vibrations decrease the contrast or 
totally ” wash-out” the interference fringes. For that rea¬ 
son good shielding and damping is required. 

However, the technological progress allows for particle 
detectors with a high spatial and temporal resolution. In 
this paper we demonstrated that the additional temporal 
information can be used to determine and decrease the 
effect of dephasing significantly even if the interferometer 
is under the influence of a multi-frequency perturbation 
and broad-band noise. We presented a theoretical ap¬ 
proach that uses second order correlation theory to ana¬ 
lyze an interference pattern disturbed by a spectrum of 
various frequencies with different amplitudes and phases. 
Thereby, the specific composition of the perturbation fre¬ 
quencies is relevant. It determines, whether the exact 
solution is needed, or an approximate solution, requiring 
less computational effort, is sufficient. 

Our method was verified experimentally by electron 
interference patterns that were artificially perturbed and 
” washed-out” by several superposed frequencies. The re¬ 
sulting correlation data for two and three frequency per¬ 
turbations and a perturbation noise spectrum matched 
well with our theoretical description. We demonstrated 
the reconstruction of the unperturbed interference pat¬ 
tern with the obtained information about the perturba¬ 
tion frequencies and amplitudes. 


The method can be applied if the perturbation frequen¬ 
cies are lower than the average count rate of the detected 
particles. This requirement is fulfilled for most interfer¬ 
ometry experiments, where a count rate of several kHz 
and perturbation frequencies of a few hundred Hz are 
common. The obtained information from our analysis 
allows for the optimization of damping and shielding in¬ 
stallations. Our technique has a potential application for 
(mobile) interferometers in a noisy environment and in 
sensor technology. The susceptibility of electron interfer¬ 
ometers may be used for the construction of a sensor with 
a high sensitivity for electromagnetic and vibrational fre¬ 
quencies. The resolution of such a device would presum¬ 
ably only be limited by the count rate of the particles 
and the distance of the interference fringes. It can be 
increased by the use of heavier particles, such as ions, or 
a larger beam path separation [34] . 

Our method identifies and analyzes multi-frequency 
noise in interferometric measurements. It is applicable 
for a broad variety of matter-wave interferometers, where 
the reduction of dephasing and the proof of the wave 
nature of particles is required in a noisy environment 
with multi-frequency perturbations. 
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